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Summary. — The Feynman-Metropolis- Teller treatment for compressed atoms [1] 
is here reconsidered in the framework of the relativistic generalised Fermi- Thomas 
model, obtained by Ruffini et al. [13]. Physical properties of a zero temperature 
plasma is thus investigated and the resulting equation of state, which keeps into 
account quantum, relativistic and electromagnetic effects, is applied to the study of 
equilibrium configurations of relativistic White Dwarfs. 

It is shown that numerical evaluation of such configuration leads, for the same 
central density p c , to smaller values of radius R and of mass M than in the classical 
works of Chandrasekhar [14] and Salpeter [10] , the deviations being most marked at 
the lowest densities (up to 30% from the Chandrasekhar model and 10% from the 
Salpeter one for p c ~ 10 6 g/cm 3 , corresponding to M ~ 0.2Mq). 
At high densities we considered the occurrence of inverse beta decays, whose effect 
is to introduce gravitational instability of the configurations. We consequently find 
the maximum mass of White Dwarfs, which, for an Oxygen and an Iron WD, is 
respectively 1.365M0 and 1.O63M0. 
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1. Introduction 

That white dwarfs cannot have an arbitrary large mass and that a critical mass 
against gravitational collapse must exist in their configuration of equilibrium reachable 
by increasing their central densities was clearly understood in the thirties independently 
by Chandrasekhar [2, 3] and by Landau. [5] 

Actually an interesting alternative approach, with respect to the evaluation of the 
interaction term between the electrons and the nuclei, was advanced two years earlier than 
the work of Chandrasekhar and Landau by the Soviet physicist Yacov Ilich Frenkel.[4] 
This paper, quoted by Chandrasekhar in his classic book, [14] has been in the past and 
is still today surprisingly ignored. Frenkel proposed to use a relativistic Thomas-Fermi 
model within a Wigner-Seitz approximation in order to describe stellar matter. The 
astrophysical motivations were not clear to Frenkel, though the theoretical formulation 
he proposed would have deserved a much more thorough examination. 
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An important turning point came in this field due to the fundamental work of Feyn- 
man, Metropolis and Teller. While at los Alamos, they considered the equation of state 
of compressed matter [1]. 

It was Ed Salpeter[9, 10] who first applied the Frenkel equations and the Feynman- 
Metropolis-Teller[l] approach to the treatment of white dwarfs: Salpeter's treatment 
was somewhat affected by the intricacies of the numerical integrations and the drastic 
approximations he adopted in the numerical solutions. 

Today with relatively simple numerical computations, we can reconsider the prob- 
lem by using a relativistic generalisation of gravitational Fermi-Thomas like equations 
[11, 13, 6] and present a comprehensive treatment by first considering the issue of the 
electromagnetic interactions between nuclei and electrons and turning later to the issue 
of inverse beta decays. The treatment presented here leads to the disappearance of all 
scaling laws with the chemical composition of the stars, which have been commonly as- 
sumed in all treatments up to now. A new value of the critical mass of white dwarfs is 
obtained, smaller than the one originally proposed by Chandrasekhar and by Salpeter, 
and one which strongly depends on the chemical composition of the star. These results 
appear in principle to be observationally testable. Particularly attractive is the possibil- 
ity of explaining the observed masses of binary pulsars: these systems are clearly neutron 
stars which could be formed by the onset of instability of stellar cores very close to the 
critical mass of white dwarfs as presented by Chandrasekhar [15]. The evaluation of their 
masses will be clearly affected by our computations. 

Here, following Fcrrcrinho et al. [11] and Ruffini et al. [13] we reconsider the Feynman- 
Mctropolis-Tcllcr[l] approach in the framework of a generalisation of the Fermi-Thomas 
model, which is here briefly derived and solved (section 2) for different atomic species in 
different states of compression. Results so obtained are applied to determine the equation 
of state (section 3) of cold degenerate stars (White Dwarfs) to show how the relativistic 
and electromagnetic effects acting on microscopic scale affects the overall structure of the 
star. Finally, the equation of state obtained is used to determine (section 4) the equilib- 
rium configurations of White Dwarfs and in particular the maximum mass configurations 
for different chemical compositions. A comparison is made with the classical works of 
Chandrasekhar and Salpeter to show the differences introduced by the use of the new 
equation of state, and to show how that they correspond respectively to a zero-order and 
first-order approximation of the complete treatment here presented. 

2. The generalised Fermi-Thomas model 

Let us consider the spherically symmetric problem of a nucleus with Z protons and 
A nucleons interacting with a fully degenerate gas of Z electrons. 

The fundamental equation of electrostatics for this problem reads, with usual meaning 
of symbols, 

(1) AV(r) = Anen e (r) - Airen p (r) 

where we introduced the number density of electrons n e (r) and of protons n p (r). 

In the equation 1 the quantities V(r) and n(r) are of course not independent, being 
related by the equilibrium condition for a relativistic gas in a coulomb external potential 
(see [16]) 



(2) 
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where the name Ep stands for Fermi- Thomas chemical potential or Fermi Energy of the 
electrons. 

To put the equation 1 in a simple and adimcnsional form we introduce the new 
function $(r), related to the coulomb potential by the formula 

(3) $(r) = V{r) + E F /e 

and the corresponding adimcnsional function x, implicitly defined by 

Ze X 



(4) 



$(r) 



Furthermore we introduce the new independent variable x, obtained by the radius r 
with the relation r = bx, where we put 



(5) 



b = (3?r) 
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Using 3 it is possible to put eq. 2 in the form 
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which, using 4, becomes 

(7) p F = 2mc 
where 
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So we obtained an expression for the Fermi momentum in terms of adimcnsional quan- 
tities, and if we remember the relation between the Fermi momentum and the number 
density of a degenerate fermion gas 



(9) 

we obtain the following expression 

Z 



n e = 
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Let us try to express the second term of the right-hand side of eq.l in terms of 
adimensional quantities: we shall assume an homogeneous spherical nucleus, with a 
radius given by the approximate formula 



(11) 
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Fig. 1. - Results of numerical integration of generalised Fermi- Thomas equation for Helium 
(3 columns on left) and Iron. For each atomic specie we give the values of the adimensional 
function y f° r different atomic sizes (i.e. for different xo&nd the corresponding adimensional 
Fermi momentum. 
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The number density of protons is therefore 

3Ze 



n v = 
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Finally we can write eq.l in the form 
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where x nuc is the adimensional size of the nucleus (r nuc = bx nuc ). 

Equation 13 is what we call ' Generalised adimensional Fermi- Thomas equation^)' . 

( 1 ) Classical Fermi-Thomas model is easily recovered in the limit (-^-) 4 ^ 3 f — > and x nuc — > 0, 
i.e. in the limit of non-relativistic expression for the momentum and of point-like nucleus. 
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The first initial condition for this equation follows from the fact that x oc r$ and 
therefore \ 0) an d so 

(14) X (0) = 

The second condition comes from the normalisation condition 



(15) 
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with r = bx atom size. Developing this formula we have 

3Z 
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which gives the simple relation 

(17) N = Z[x aX '(x )-x(xo) + l} 

For a neutral atom we have N = Z and the condition 17 reads simply 

(18) xox'(xo) = x{x ) 

Note that the physical quantities of interest, such as the coulomb potential and the 
density of electrons do not show any singularity in the center, neither on the border 
of the nucleus, being dependent just on the function X an d his first derivative. At the 
opposite the only discontinuity appears in the second derivative of x due to our rough 
assumption of homogeneous spherical nucleus. 

We show in figure 2 a numerical integration of 13 for a |g-Fe atom, in several states 
of compression. 



3. The equation of state 

The study of the previous section allows us to determine, under suitable hypothesis, 
the equation of state of compressed matter, by simply computing for each value of the 
atomic compression parameter xq the corresponding values of pressure and mass density. 
The origin of these two physical quantities is evidently different, the first being generated 
from the degenerate gas of electrons and the second being essentially given by the mass 
density of nuclei. We shall assume that nuclei arc arranged in a Wigner-Seitz lattice, 
each cell being filled by a relativistic gas of degenerate electrons. As shown by Salpeter 
[9] the shape of the lattice cell is in first approximation unimportant, so we'll assume a 
spherical cell, with a radius equal to the atomic radius. The average density of cells is 
thus equal to 

Am 3Am 

(19) p(x ) = 



Vcell 4tt(&X ) 3 

where A and m are respectively the number and the mass of nucleons. 
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Fig. 2. - Solution of the generalised Fermi- Thomas adimensional equation for a Fe atom in 
several states of compression. 



In this scheme also the pressure is easily determined, being simply the pressure gen- 
erated by a Fermi gas of number density equal to the number density of electrons at the 
border of the atom n e (x ), which can be expressed in terms of adimensional quantities 
as follows 



(20) nM = ^{—) [ 1+ UJ — _ 

In tab. 2 we list the resulting equation of state for different atomic species. 

We would like to stress the deep difference between our model and the Chandrasekhar's 
one, evidently corresponding to an electron gas uniformly distributed within the cell, by 
showing the behaviour (see figure 4) of the ratio R(xo) between the density at the border 
and the average number density, given by 
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Fig. 3. - Equation of state obtained for Iron (the two columns on the right) and for Helium. 
Pressure and mass densities can be obtained directly the quantities in tab.l (see text). 



(21) R = ^pf = i (xox(-o)) 3 

n{x ) 3 

The analysis of this plot clarifies some important features: 

• For high values of xq, i.e. for low densities the ratio R reduces to the non-relativistic 
(dashed line) expression 



( Zcr ^ 
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(22) Rnonrel(xo) = ^ (x x(a;o)) 3/2 



• For lower value of xo , corresponding to values of density up to an important fraction 
of nuclear density one gets the ultra-relativistic value 



(23) Rultrarel = ^ {^~^ *( X o) 3 

It can be shown that scaling laws apply in this case and that the ratio is in effect 
independent from the compression parameter (see [12]). 

• Finally for x = x nuc we recover the uniformity of the electron distribution. 

But our treatment allows us not only to remove the zero-order approximation of 
Chandrasckhar (uniformly distributed electrons) but also the first-order approximation 
of Salpeter who could not solve the TF equation for each value of Z and £0, as is done 
here, and introduced a development of the number density of electrons 

(24) n(x) = n + e(x) 

which is evidently only allowed when the electrostatic interaction within nuclei and elec- 
trons is small respect to the Fermi Energy of electrons, i.e. in the limit of high densities. 

4. Equilibrium configurations of white dwarfs 

We are now able to compute equilibrium configurations of white dwarfs. It is clear 
that we will obtain, starting with the same central density, less massive configurations 
than those obtained by both Chandrasekhar and Salpeter. That is due to the fact that 
for a fixed mass density i.e. for a fixed value of A and of the compression parameter, our 
pressure is systematically lower being generated by a shell of electrons whose density is, 
as we shown, always lower than the average one and eventually very far from the mean 
value. 

Numerical integration of the models are here performed for stars of Helium and Iron 
in the framework of a Newtonian theory of gravitation and General Relativity. In the 
last case we have the usual TOV equation of hydrostatic equilibrium 

(25a) = 4nr 2 e(r) 

dr 

( GM(r) 

(25b ) * - JziAL^l 

K ' dr C 2 1 2GM{r) 

c' 2 r 

and we add to these equations the numerical equation of state obtained in the previous 
section. 

Results of integrations, plotted in figg.5 and 6, show that deviations from classical 
models become very important at the lower densities. 
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Fig. 4. - Ratio R(xo) between density at the border and average density for Fe, O and He atoms. 
It is shown also the corresponding results obtained with a non-relativistic treatment. 



The mass for an Iron star corresponding to a central density of the order of 10 6 <?/cm 3 
is about 30% smaller than what predicted from the Chandrasekhar model and about 10% 
smaller than the Salpeter one. Of course deviations are most marked for atomic specie 
with high Z, being related, as seen, to Coulomb corrections to the equation of state. 
Numerical integration of equilibrium configurations can be computed up to arbitrarily 
high central densities, but to make the model physically reasonable we have to make 
allowance to the occurrence of nuclear reactions. In particular when the Fermi Energy of 
electrons is sufficiently high they can undergo inverse beta decay, consisting in a reaction 
of the type 

e~ + p — > n + v e 

Using the recent (1993) nuclear data in Audi et al [17] we can find for each element 
the beta decays energies, and thus the critical density at which the element considered 
undergoes inverse beta decay. As first shown in Harrison et al [18], the introduction of 
nuclear reactions leads to the existence of a maximum mass which separate stable from 
unstable configurations (see [19], based on the fundamental works of Chandrasekhar 
[20] [21]). Using this data, we stop the curve when the central density is equal to the 
critical density. The configuration wich correspond to the last point is the maximum 
mass configuration. For an Oxygen and an Iron WD, we find respectively 1.365M Q and 
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Fig. 5. - Equilibrium configurations curve for Oxygen WD in the M/Mq — p c plane, obtained 
in Newtonian theory and General relativity. For comparison we show the results obtained with 
Chandrasekhar model. 



1.063M Q . These values differ roughly about 5% from the Salpeter configurations and 
10% from those evaluated by Chandasekhar (see figures). 

5. Conclusions 

The derivation of the relativistic generalised Fermi-Thomas model is here presented, 
and a FMT treatment of the compressed atom is applied. Numerical examples are 
given of compressed atomic configurations and an equation of state is derived for a 
zero temperature plasma, removing some unneeded approximations applied in classical 
papers. Results are applied to the problem of equilibrium configurations of White Dwarfs, 
described as cold degenerate stars and a class of configurations is obtained (for the same 
central densities) less massive than those predicted in the models of Chandrasekhar and 
Salpeter; in particular we give an estimate of the maximum mass for Oxygen and Iron 
WD. 
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